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Abstract 

The present paper studies a large class of temperature-dependent 
probability distributions and shows that entropy and energy can be 
defined in such a way that these probability distributions are the equi- 
librium states of a generalized thermostatistics, which is obtained from 
the standard formalism by deformation of exponential and logarith- 
mic functions. Since this procedure is non-unique, specific choices 
are motivated by showing that the resulting theory is well-behaved. 
In particular, with the choices made in the present paper, the equilib- 
rium state of any system with a finite number of degrees of freedom is, 
automatically, thermodynamically stable and satisfies the variational 
principle. 

The equilibrium probability distribution of open systems deviates 
generically from the Boltzmann-Gibbs distribution. If the interac- 
tion with the environment is not too strong then one can expect that 
a slight deformation of the exponential function, appearing in the 
Boltzmann-Gibbs distribution, can reproduce the observed distribu- 
tion with its temperature dependence. An example of a system, where 
this statement holds, is a single spin of the Ising chain. However, be- 
cause all systems of the present generalized thermostatistics are au- 
tomatically stable, one must not expect that all open systems can be 
described in this way. Indeed, systems exhibiting a phase transition in 
their thermodynamic limit, can be unstable even when the interaction 
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with the environment is weak. Therefore, their equilibrium probabil- 
ity distribution cannot be described by a simple deformation of the 
Boltzmann-Gibbs distribution. In order to be able to handle such sys- 
tems as well the second part of the paper discusses a further extension 
of the class of probability distributions using mean- field techniques. 

Connections are discussed that exist between the present formal- 
ism, Tsallis' thermostatistics, and the superstatistics formalism of 
Beck and Cohen. In particular, the present generalization sheds some 
light onto the historical development of the Tsallis formalism. It is 
pointed out that temperature dependence of the Tsallis distribution 
has hardly been verified experimentally. 

Keywords Generalized thermostatistics, mean-field theory, de- 
formed logarithmic and exponential functions, Ising chain, non-exten- 
sive thermostatistics, superstatistics 

1 Introduction 

Statistical physics is primarily concerned with (nearly) closed systems. A cor- 
rect description of open systems, i.e. systems weakly interacting with their 
environment, requires that the latter is taken into account explicitly. The 
present paper attempts to generalize the thermostatistical formalism to open 
systems without the need for including the environment. Such a generaliza- 
tion may fill the gap between the narrow context for which the theory of sta- 
tistical physics is developed and the vast domain of its applications. Indeed, 
often used tools, like temperature-dependent Hamiltonians and mean-field 
approximations, are not so well founded in the standard formalism. The 
main part of the paper explores how far one can go in modeling the tem- 
perature dependence of open systems. In doing so, Hamiltonians are not 
allowed to depend on temperature. Rather, the Boltzmann-Gibbs distribu- 
tion is replaced by more general probability distributions and the microscopic 
definition of entropy is modified accordingly. Only the last sections include 
temperature-dependent Hamiltonians as a further extension of the formalism, 
more suited for unstable systems. 

A general formalism like the one developed here helps also to improve 
understanding of more specific generalizations of the Boltzmann-Gibbs for- 
malism. The present paper discusses in particular the formalism of non- 
extensive thermostatistics, first introduced in pQ. By looking from an eagle's 
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point of view several peculiarities of this formalism get a new interpreta- 
tion. Let me mention the role of escort probabilities, historical problems 
with average energy without proper normalization and with the definition of 
temperature, and the occurrence of thermodynamic instabilities. Also the 
relation with the recently introduced superstatistics j2] is discussed. In this 
way the present paper contributes to clarify the relation between mainstream 
statistical mechanics and its generalizations. 

It is of course very important that the present generalized formalism has 
physical applications. Nevertheless, no such applications will be discussed. 
The literature of non-extensive thermostatistics claims a large number of 
applications — see e.g. |3j. It is argued below that these applications are 
candidate applications of the present formalism. However, in most, if not all, 
of these papers the authors verify the shape of the equilibrium probability 
distribution, but not its temperature dependence. The present paper stresses 
the point that in a predictive theory of thermostatistics both aspects are 
equally important. Let us now explain this point. 

Gibbs' postulate implies that microstates with higher energy have a lower 
probability, according to the formula 



In this expression pk is the probability of the microstate labeled k. It sat- 
isfies pk > and J2kP k = 1- -fffc is the corresponding energy and T is the 
temperature (Boltzmann's constant is taken equal to 1). 

Mathematically, it is clear that any probability distribution function (pdf ) 
with non- vanishing probabilities can be written in the form (JTJ. This requires 
only an appropriate definition of the notion of energy. Indeed, given p^, one 
can put 



But, the Boltzmann-Gibbs distribution contains more information than 
the mere statement that high energy states have an exponentially small prob- 
ability. It also predicts how probabilities pk change with temperature T. 
Once temperature is taken into account a definition of energy, like (|5]L is not 
acceptable because in general energy levels Hk, defined in such a way, are 
temperature dependent. One concludes that the Boltzmann-Gibbs distribu- 
tion (JTJ) is very specific once temperature dependence is taken into account. 
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In the terminology of the present paper it is a temperature-dependent prob- 
ability distribution. 

The main tool of the paper is the notion of K-deformed logarithmic and 
exponential functions, a definition of which is given in Section III. The idea 
of considering density functions as deformed exponential functions goes back 
to jl] . The name of K-deformed functions has been used explicitly in [5. . One 
could of course replace the Boltzmann-Gibbs distribution (JTJ) by an arbitrary 
temperature-dependent probability distribution. However, chances are small 
that in such a generality one could develop a complete thermostatistical the- 
ory. Therefore it is obvious to require that the function, which replaces the 
exponential function in JI}, satisfies some additional properties. In the first 
place it should be an increasing function because states with a higher energy 
should be less probable. An additional property of the exponential function 
is convexity (i.e. the second derivative is positive). Functions satisfying both 
properties (plus some technical conditions) are called K-deformed exponen- 
tials and are denoted exp K (x). Doing so has the advantage that well-known 
expressions from standard statistical physics are recognized on sight in the 
generalized formalism. 

The next section starts with a discussion about the relevance of pdfs, 
other than that of Boltzmann-Gibbs. In the third section a new approach, 
based on K-deformed logarithmic and exponential functions, is proposed. 
Section 4 discusses the Ising chain as an example motivating the generalized 
formalism following in section 5. Further examples are given in section 6. 
This includes a discussion of Tsallis' thermostatistics. Generalized mean-field 
models are treated in section 8, using the notion of probability dependent 
variables introduced in section 7. Finally, a short discussion follows in section 
9. The paper ends with two appendices, one about the Ising model, the other 
about thermodynamic stability. 

2 Evidence of limited validity of the Boltz- 
mann-Gibbs distribution 

Consider an open system described by probabilities p^, where k labels the 
microstates of the system, and 7 the microstates of the environment. Assume 
these probabilities are of the form (JTJ, i.e. there exist energy levels such 
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that 

Pkl = 'zjrj ex P( _i W T )- ( 3 ) 
Write these energies in the form 

H kl = H k + h; + V kT (4) 

Then the reduced system without the environment is described by the prob- 
abilities 

p k = $> fc7 = Yjf) eM ~ {Hk + (Vfc) * )/T) ' (5) 

7 

with the non-linear Kolmogorov-Nagumo averages 0[7| (Vk)* given by 

(Vk)* = -Tin ^ P ;exp(-y fc7 /T)j , (6) 

with 

P^y^rfM-^/n (7) 

and with appropriate normalizations Y(T) and Y*(T). Clearly, the reduced 
probabilities p k are not of the form (JTJ), since the effective energies H k + (Vk)* 
depend on temperature. The only exception is when the interaction V kl be- 
tween system and environment vanishes. One concludes that the Boltzmann- 
Gibbs distribution (JTJ can only be universally correct in the limit that in- 
teractions of the system with its environment are negligible. Note that for 
systems of statistical mechanics the interactions cannot vanish because then 
the system is no longer open. Indeed, the dynamics of closed systems predicts 
conservation of energy. As a consequence, the state of the isolated system 
is not described by (JIJ) but by a probability distribution concentrating on a 
single value of the energy. 

The above reasoning shows that temperature-dependent effective Hamil- 
tonians are rather generic. They occur in many situations and are used so 
frequently in applied statistical physics that it is a labor of Sisyphus to quote 
all occurrences. Nevertheless, they do not receive much attention in the 
theoretical foundation of statistical physics. 
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In the above discussion of open systems the assumption is made that the 
totality of system plus environment is correctly described by the Boltzmann- 
Gibbs distribution. This is not a severe limitation. Indeed, a characteristic 
of systems of statistical physics with many degrees of freedom is that the 
number of energy levels Hk in an interval [E, E + AE] increases exponentially 
with increasing energy E. As a consequence, the typical equilibrium pdf p 
has the property that the only microstates that contribute significantly have 
energy levels Ek approximately equal to the average value U = J2kPkHk, 
in the sense that the sum of pu over these states is close to 1. This is the 
basis for the equivalence of ensembles — see jHJ El EH] for a mathematical 
formulation of this property. In information theory ^T] this property is 
known as equipartition theorem. 

The equivalence of ensembles, when valid, implies that the actual form 
of the pdf is not very crucial. It is only needed for microstates with energy 
level Hk in the vicinity of the average energy U. In particular, most systems 
of statistical mechanics cannot be used to test the validity of the Boltzmann- 
Gibbs distribution. 

3 A new approach 

The present study starts from a family of temperature-dependent probabili- 
ties, such as could be observed for an open system. If this family is of a certain 
form then one can construct a thermostatistic formalism that produces these 
probabilities as the temperature-dependent equilibrium distribution. 

In order to generalize (0), the recently jT2] introduced notions of re- 
deformed exponential and logarithmic functions are used. A re-deformed 
logarithmic function is denoted ln K (x) and is a strictly increasing concave 
function defined for all positive x, normalized such that ln K (l) = 0, and with 
a possible divergence at x — 0, mild enough so that the integral 



converges for x — 0. This definition is not very restrictive. A simple example 
of a re-deformed logarithm is ln K (x) = Sa/x — 3. Examples of families of 
re-deformed logarithms follow below — see (J23HZI- 

The inverse function of a re-deformed logarithm is a re-deformed expo- 
nential function and is denoted exp K (a;). If x is not in the image of the 
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K-deformed logarithm then expire) is taken equal to when x is too small, 
and equal to +00 when x is too large. 

Fix now a K-deformed exponential function exp K (x). Note that through- 
out the paper this function does not depend on temperature. The obvious 
generalization of (JTJ) is now 

Pk (T)=ex PK (G(T)-H k /T), (9) 

where the function G(T) is chosen such that normalization Y2 k P k — 1 is 
satisfied. The alternative of putting normalization into the partition sum 
Z(T), i.e. 

p k (T) = ^ f - ) ex Pfi (-H k /T) (10) 

is discussed later on. Expressions (jUJ) and (fTUJ) are of course equivalent in case 
the deformed exponential is replaced by the standard exponential function. 
In fact, the difference between the two expressions is only relevant when 
temperature dependence is considered. The choice (jOJ is the simpler one 
because it yields immediately an expression for the quantity ln K (p k ), needed 
in the context of the definition of entropy (see below). 



4 Ising example 

The following example illustrates well the point of view of the present paper. 
The Ising chain in an external field can be solved exactly in the thermody- 
namic limit. The result for the occupation probability of the two levels of 
the spin at the origin is 

P±i T) = i ( 1 ± - 7 J£L=) , 

with 

y(T) = sinh (^) e? J . (12) 

In the latter expression A > is the energy level splitting due to the external 
field, J > is the interaction energy, and (3 = 1/T is the inverse temperature. 
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Figure 1: Deformed logarithm for the Ising model with J/A = 0.25, accord- 
ing to (JSTJ E21 EH}- The dotted line is the natural logarithm. 

The question under study is whether there exists a thermodynamical 
formalism which reproduces the temperature-dependent probabilities p±(T) 
without reference to the environment, which in this example is formed by all 
other spins of the Ising chain. The obvious definition of internal energy U is 



If there is no interaction with the environment (i.e. J = 0) then the pdfs 
(|11|) are of the Boltzmann-Gibbs form. In this case, Shannon's measure of 
information 



is an adequate expression for the entropy S. Indeed, it is well-known that the 
Boltzmann-Gibbs distribution is thermodynamically stable in the standard 
formalism of statistical physics based on (|14|1. For further use let us verify 
that temperature T corresponds with the thermodynamic notion of temper- 
ature. A short calculation using (jl4j) and the Boltzmann-Gibbs distribution 
gives 



U=--\p+{T)-p-{T)]A. 



(13) 




(14) 




(15) 
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so that 




(16) 



Hence, the thermodynamic relation 



is fulfilled. 

If the spin interacts with its environment (i.e. J ^ 0) then the pdfs 
(I11J) are not of the Boltzmann-Gibbs form, but can be written in the form 
at least when J/ A is not too large — see Appendix A. For the proof 
of this statement it is essential that inverse temperature can be written as 
a monotonic function of the probabilities p + and p_. Fig. 1 compares the 
deformed logarithm, used in this case, with the natural logarithm. In what 
follows will be shown that, with an obvious generalization of the definition 
of entropy ()14|). the thermodynamic relation ()17|) still holds. The definition 
of internal energy (fT3|) is not modified. 

5 Generalized thermostatistics 

Given a temperature-dependent family of probabilities Pk{T) of the form 0, 
a generalized thermostatistics, which for all temperatures yields @U) as the 
equilibrium distribution, is formulated as follows. 

Let energy as a function of temperature be defined by 




(18) 



k 



and let entropy as a function of temperature be defined by 



S{T) 




((-O)p(T) With 



(19) 



k 
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In this expression, F K (0) is a constant defined by (jSJ)- The extended class of 
entropy functionals of the form (|19|) has been studied in fT2\ IT%j . A short 
calculation gives 



§ = £[-Mp,)-f k( o)]^ 

k 

= £[-G(T) + tf fc /T-F K 
k 

1 dU 



dpk 
dT 



TdT 

to see the latter, use that 



^ dT 

k 

This proves the thermodynamic relation ()17j) . 

Next, introduce so-called escort probabilities jEE], defined by 

P k (T) = Z(T)- 1 exp> K (G(T)-H k /T) 

with exp' K (x) the derivative of exp K (x) and with 

Z(T) = J2^P' K (G(T)-H l /T) 



(20) 



(21) 



(22) 



(23) 



These escort probabilities have been introduced in jTHj for the special case of 
Tsallis' thermostatistics. From Y2 k Pk = 1 f°U° ws 



dG 
dT 



-T- 2 Y,PkiT)H k . 



(24) 



dT 



Using this result one calculates 
dU 

k 

= Z(T)T- 2 J2 p k(T) 

k 

> 0. 



-1 2 



H k -J2 p i( T ) H i 



(25) 
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Hence, energy is an increasing function of temperature, as it should be. To- 
gether with (|17p. this shows that entropy S(U) as a function of energy U is 
concave. In other words, the system is thermally stable. 

There exists also a variational principle, satisfied by the probabilities p k . 
Given any pdf r, consider the problem of maximizing information content 
{{I)) r under the constraint that average energy J2 k r k H k equals a given value 
U. Introduce Lagrange multipliers a and /3 and minimize the expression 

/3£Y fc ff fc -((J)) r -ar£V fc . (26) 

k k 

Variation w.r.t. r k gives the condition 

= PH k + \n K (r k )+F K (0)-a. (27) 

Comparison with Q shows that r = p, (3 = 1/T, and a = F K (0) + G(T), is 
a solution of this variational principle. 

Assume now that M is any other other observable of the system, e.g., 
total magnetization of a spin system. Then thermodynamic stability ^7] 
requires not only that entropy S is concave as a function of internal energy 
U but also as a function of the average value (M). That this property holds 
is shown in Appendix B. 



6 Examples 

The example of a single spin of the Ising chain has been discussed above. The 
pdfs flllj) are of the form @. Therefore, the formalism developed above is 
applicable. One concludes that the exact equilibrium distribution of a weakly 
interacting Ising chain in the thermodynamic limit, when restricted to a single 
spin variable, is an equilibrium distribution of generalized thermostatistics. 

6.1 Tsallis' thermostatistics 

Make the following specific choice of K-deformed logarithm 

ln K (x) = (x 9 " 1 - 1) , (28) 

o — l 
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where q is a parameter, which for technical reasons should lie between and 
2. The inverse function is 



exp re (xj 



1 ' 

1 H x 

Q 



1/(9-1) 



(29) 



The notation [x} + = max{0, x} is used. A short calculation shows that (|T9"|) 
reduces to 



(30) 



which is Tsallis' g-entropy PQ. The temperature-dependent probabilities p k , 
given by equation (jUJ), read 



1 + £Z1 ( G (T) - * 



1 1/(9-1) 



(31) 



The distribution introduced by Tsallis in its original version, reads in 
the present notations 



Pk = -[l-(3*{q-l)H k 



1/(9-1) 



(32) 



with Z a normalization constant, and (3* a Lagrange multiplier related to 
inverse temperature. It is of course possible to write (|31|) in the form of (|32|). 
The identification gives 



(3* 
Z 



Tq+(q-l)G(T) 

! i i/d-9) 
1 + G(T) 



(33) 



Note that G(T) can be eliminated to obtain 



1 



ft* = —z g -\ 

qT 



(34) 



Clearly, the Lagrange multiplier (5* does not have its usual value of 1/T, but 
depends in a rather complicated way on temperature. This difficulty was 
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not understood in the early days of Tsallis' thermostatistics. In fact, it took 
about ten years [TH] before the question of defining temperature was settled. 

It is interesting to sketch briefly the evolution that the formalism has 
undergone. Soon after the original paper [T] appeared, it was proposed [TH] 
to replace the constraint J2k rk ^ k = ^> usec ^ ^ n the variational principle, 
with a constraint of the form 

J24H k = U. (35) 

k 

The effect of this change is a pdf of the form with q — 1 replaced by 
1 — q. The main advantage of the above constraint is that the usual relation 
(3 = 1/T between Lagrange multiplier and thermodynamic temperature T, 
defined by (|17j) is restored, however, at the cost of changing the definition of 
thermodynamic energy U. Both papers jU EI] have been criticized E] > 
the first because of the apparent lack of thermodynamic stability, the latter 
because internal energy U, as given by (|3*5j). is not the average of the energy 
levels Hk- Later attempts to cure the formalism [THIHE] (introducing escort 
probabilities and changing the definition of temperature) solved the problems 
only partly. For a recent discussion, see |22j . 

6.2 Superstatistics 

In a recent proposal [2S1 EH 121 E2J the Boltzmann-Gibbs formalism is gen- 
eralized by replacing the Boltzmann factor exp(— f}E) by an average over a 
range of inverse temperatures f3 

POO 

B(E)= / d(3f((3)e^ E . (36) 
Jo 

The resulting formalism is called superstatistics. With an appropriate choice 
of distribution /(/?) the Tsallis distribution (|3*2|) is recovered as 

Vk = ^B(H k ) 
1 f°° 

= zJo d(3f ^ e (37) 

with Z = J2k B(H k ) (adapting notations to those of the present paper). The 
argument in favor of ()36| l37j) is that in non-equilibrium systems thermo- 
dynamic variables like temperature may fluctuate in space and time. This 
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argument is not so different from the open systems argument, given above, 
that due to interactions with the environment even the equilibrium distribu- 
tion will deviate from Boltzmann-Gibbs. 

Note that all probability distributions of superstatistics can be written as 
distributions of the present generalized thermostatistics. In other words, any 
distribution of the form (|37|)) can be written as (JSJ) with an appropriate choice 
of K-deformed exponential function. Indeed, fix some positive temperature 
T and let 

poo 

exp K (x)= / d(3g((3)e T ^ (38) 
Jo 

where g(/3) is a probability distribution still to be determined. The function 
exp K (x) is increasing and convex. It satisfies exp K (0) = 1. One has 

poo poo poo 

/ dxexp K (-x) = / d(3g{(3) / dxe~ Tfix 
Jo Jo Jo 

= [ ipm w (39) 

The latter integral is convergent if the distribution g(0) stays away from 
(3 = 0. Then exp K (x) satisfies all requirements to be a deformed exponential 
function in the sense of ^2j. Now use (jHJ) to calculate 

POO 

p k = / d/W)e T/3G(T) -^. (40) 
Jo 

This expression coincides with (|3*7jl provided 

9ifi) = ^f(P)e~ XP (41) 
with A the unique solution of 

d(3f((3)e- x P = Z. (42) 



The normalization function G{T) is then given by G(T) = X/T. 

The above arguments show that any probability distribution produced 
by the formalism of superstatistics can also be produced by the present for- 
malism of generalized thermostatistics. Is it possible to distinguish the two 
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formalisms by applying them to problems of physics? The obvious answer 
is that one should look to the temperature dependence of the equilibrium 
probability distributions and compare them with experimental data. How- 
ever, superstatistics is meant to be valid for systems not in equilibrium, with 
a fluctuating temperature. It is not at all clear how the distribution of inverse 
temperatures f(/3) should depend on some average temperature. This point 
deserves further investigation. 

The inverse function of exp re (;c) given by (|3*8)l is denoted ln K (x), as before. 
It can be used to calculate an entropy functional ((I)) p by means of (fT^j). 
Up to a prefactor, this is exactly the entropy functional introduced in |2o] . 
in the first of the two cases considered there. Indeed, the expression, found 
there, reads 

rPk 

S(P) = J2 dy(a + E(y)) (43) 
k Jo 

with a a constant and E(y) the inverse of the function Z~ 1 B(E). In the 
present notations is 

ii?(£)=exp K (i(A-£)). (44) 

Hence the inverse function is given by 

E(y) = \-T\n K (y). (45) 

One concludes that 

fPk 

S(p) = J2 dy(a + \-T\n K (y)). (46) 
k Jo 

This equals T times ({I)) p provided that a is chosen equal to —A — TF K (0). 

Note that the formalism developed in 26J deviates from the one presented 
here because energy constraints other than the standard (|18|) are considered. 
Such constraints are considered in the final part of this paper. 

6.3 Kaniadakis' entropy functional 

Another example of /t-deformed logarithm has been introduced by Kani- 
adakis jSJ I2ZI , namely 

Ms) = ^.(^- » ( 4 7) 



15 



with k a parameter with value between -1 and 1. The corresponding deformed 
exponential function is 



exp K (x) = 



KX + Vl-f K 2 X 2 



(48) 



With these definitions the entropy functional (fTHj) becomes 



The corresponding equilibrium distribution is 



Pk = 



(«[G(T) - tf fc /T] + v/l + K 2[G(T)-/f fc /r]2) 



(50) 



This distribution differs from the one studied in [3] because here normaliza- 
tion is not a prefactor. 

7 Probability-dependent variables 

In this final part of the paper the class of temperature dependent probabilities 
(jUJ) is enlarged. One cannot expect that all open systems, weakly coupled 
with their environment, are described by equilibrium distributions of the 
form 0. Indeed, consider e.g. an open system consisting of a finite part of a 
2-dimensional Ising lattice with nearest-neighbor interaction constant J > 0, 
in absence of an external field. Independent of how small J is there is always 
a critical temperature T c (proportional to J) where the Ising lattice becomes 
unstable. This will be reflected by an instability of the entropy S(U, B) of 
the open system. The idea is now to capture the essence of such instabilities 
by elements coming from mean-field theory. 

Temperature-dependent variables are replaced by probability-dependent 
variables, which is a new notion that is introduced here. By definition, 
a probability-dependent variable is a function fk{x) depending on the mi- 
crostate k and on the probability x, which varies between and 1. Its 
average or expectation is defined by 




(51) 
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Note that this expression has the essential characteristics of a statistical 
average. It is linear, positive on positive functions, and normalized to one. 

The most obvious example of a probability-dependent variable is infor- 
mation content, also called entropy functional. Indeed, Shannon's measure 
of information I(jp) can be written as 



Hp) = -^PkMPk) 

k 

= ~z2 dx ( ln ( x ) + 1 

u Jo 



= «/>*>„ (52) 

with the functions Ik{x) defined by 

Ik(x) = — ln(x) — 1, independent of k. (53) 

In fact, this is the reason why f)19j) was written in the form {{I)) p , with the 
probability-dependent variable / given by 

Ik( x ) — ~ hi K (x) — F K (0), independent of k. (54) 

Notice that also energy may be a probability-dependent variable. Indeed, 
consider the following well-known mean-field Hamiltonian 

H = -Ja(a) p — ha. (55) 

In this expression J > and h are model parameters, a is a spin variable 
taking on the value ±1 with probability p±. The average value of the Hamil- 
tonian H can be written as 

(H) P = -J(a) 2 p -h(a) p 

= m) P , (se) 

with potentials 4>±{x) given by 

(j) ± {x) = -AJx + J^h. (57) 

This shows that the energy of the mean-field model (jo3|) is a probability- 
dependent variable. 
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A straightforward generalization of (J57|) is to consider potentials of the 
form 

0fcO) = -JkQX 11 ' 1 - h k , (58) 

where q is a fixed number, and Jk > and hk are model parameters. The 
corresponding average energy equals 

m) P = -J2 J *pt-J2 h ><PK- ( 59 ) 

Note that the non-linear constraint (j35|) coincides with this expression, with 
Jk < and hk = 0. Hence the critique of |20] . that (|35|) lacks normaliza- 
tion, does not apply. Indeed, in the present context the average {(<f>)) p is 
normalized to 1. 



8 Generalized mean-field theory 

The class of temperature-dependent pdfs JUJ) can now be enlarged by allowing 
all distributions of the form 

Pk = exp K (G(T)-MPk)/T), (60) 

where ipk{ x ) is a probability-dependent potential. The question, discussed be- 
low, is whether there exists a generalized mean-field model with probability- 
dependent potentials 4>k{x), for which the equilibrium distributions are given 
by (|60p. The relation between the potentials <pk(%) an d ipk{ x ) is allowed to 
be non-trivial. 

Let entropy S be defined as before by S = ((I)) p - It is immediately clear 
that internal energy U, defined by U — ((ip)) p , still satisfies (fT7j) . However, 
the stability requirement that energy U is an increasing function of T is not 
always fulfilled. One can show, e.g., that it holds when the high temperature 
condition 

T > — Vfe(P*) exp' K (G(T) - MPk)/T) (61) 

is satisfied for all k. It is not a surprise that in these mean-field models a 
thermodynamic instability can occur at low temperatures. 
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The variational principle is straightforward. Consider the problem of 
maximizing ((I)) r under the constraint that ((ip))r — U. This leads to the 
Lagrange problem of minimizing 

k 

Variation w.r.t. r k leads to the equation 

P^pk(r k ) - I k (r k ) = a. (63) 

This equation coincides with (|50|l. when r = p, (3 = 1/T, and ot(T) — F K (0) = 
G(T). 

However, the above variational calculation is not what one does tradition- 
ally with the mean-field model (|55j). Rather, one first studies the Hamiltonian 
if = — Jmer — /ur, where m is a constant which is taken equal to (a) p only 
after finishing the calculation of thermodynamic equilibrium. In order to 
generalize this procedure, consider the problem of optimizing ((I)) r under 
the constraint 

Y,nE k = u, (64) 

k 

with constant energy levels E k . At the end they are taken equal to ip k {p k ). 
The corresponding Lagrange problem leads to the equation 

f3E k - I k {r k ) = a. (65) 

It is clear that p k given by (JBTIj) is the solution of this problem. However, 
(J65|) is not the equation corresponding with the problem of maximizing ((I)) r 
given the constraint ((if))) r = U, but rather that corresponding with a con- 
straint of the form (((f))) r = U*, where the 0-potentials satisfy 

1 f y 

^kiy) = - / (pk(x)dx. (66) 

y Jo 

In the specific case of the mean-field Hamiltonian ([55)1 . in combination with 
Shannon's measure of information, ()65|) reduces to the well known mean-field 
equation 

m = tanh {(3{Jm + h)) with m = r + — r_. (67) 
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The expression 

(3j2 r k E k-((I))r (68) 

k 

should be convex, with minimum at r = p. The appropriate tool to investi- 
gate this point is relative information content, also called relative entropy or 
Kullback-Leibler distance. Its generalized definition is (see |14j ) 

I{r\\p) = I(r)-I(p)+Y J (r k -Pk)h(Pk)- (69) 

k 

Using (J33J) this becomes 

dx (\n K (x) - ln«(p fc )) . (70) 

A short calculation gives 

P^kEu - ((I)) r = (3Y,PkEk - «/»p + I(r\\p). (71) 

k k 

Because J(r||p) is a convex function, which is minimal if and only if r = p, 
this shows that, if a pdf p exists, solving ([60)1. then it minimizes ()68|) . This 
does not mean that the pdf p is the only solution of (J60|) . Indeed, we know 
that in the mean-field model (J55|) there can exist three solutions. But the 
corresponding energy levels Ek are different, so that each time (|68|) is minimal 
at r = p. 

Let us shortly return to the Tsallis formalism. Expression (|H2*|) gives 
the temperature dependence of the pdf as found in the context of Tsallis' 
thermostatistics. It is straightforward to write these in the form (|60|). Indeed, 
let ipk{x) be given by 

M$ = Hkq^ 1 (72) 

(note that this is (J58|) with some change of notation) and let \n K (x) be given 
by (EEI- Then reduces to with G(T) given by 

G ^ = ^TT ~ l ) » ( 73 ) 

and (3* = 1/T. This means that Tsallis' non-extensive thermostatistics is a 
special case of the present generalized mean- field theory. 
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9 Discussion 



This paper proposes a generalized thermostatistics obtained by replacing ex- 
ponential and logarithmic functions in the Boltzmann-Gibbs distribution, 
respectively Shannon's entropy functional, by /t-deformed functions. The 
latter are rather arbitrary functions satisfying a minimal number of require- 
ments. In the resulting formalism thermodynamic stability is automatically 
satisfied. 

New in the present approach is the emphasis on the temperature depen- 
dence of equilibrium probability distributions. Such temperature dependence 
may be obtained by calculation of an open system in the standard formalism, 
or may, in principle, originate from an experiment. Next, a thermodynamic 
formalism is sought for in which these are the equilibrium distributions. This 
way of working eliminates the need for postulating specific entropy functions 
or energy constraints, as is the typical way of doing in non-extensive ther- 
mostatistics. 

The thermodynamic relation between entropy S and internal energy U 
can be derived from the underlying microscopic statistical theory under the 
assumption that the interaction of the system with its environment is negli- 
gible. Goal of this paper is to show that this thermodynamic relation is also 
valid for systems weakly interacting with their environment. This has been 
shown for at least one example, namely a single spin interacting with its en- 
vironment. However, one cannot expect this result to be generic because the 
system together with its environment may exhibit a phase transition. In such 
a case one expects that the system might be destabilized by its environment. 

In order to be able to cope with the latter possibility a mean-field ap- 
proach has been elaborated in the second part of the paper. The concept of 
probability-dependent variables has been introduced. The result is a rigor- 
ously defined mean-field theory. Its stability properties have been analyzed 
using a generalized notion of relative entropy. Although equilibrium states 
are not necessarily unique they still satisfy a variational principle. But the 
effective energy levels involved in this variational principle depend on the 
choice of equilibrium state. 

A secondary goal of the present paper is to improve understanding of 
Tsallis' formalism of non-extensive thermostatistics. With a specific choice of 
deformed logarithmic function the entropy functional becomes that of Tsallis. 
The resulting equilibrium distribution coincides with Tsallis' distribution. If 
temperature dependence of the distributions is taken into account then the 
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Tsallis formalism describes mean-field models of the present approach. As a 
consequence, all applications of the Tsallis formalism are also applications of 
the present generalized thermostatistics. The same statement holds for the 
recently introduced formalism of superstatistics. Also the entropy functional 
proposed by Kaniadakis has been considered. The present formalism leads 
to an equilibrium distribution which differs from the one proposed by Ka- 
niadakis because here normalization is not a prefactor but enters inside the 
deformed exponential function. 

For simplicity only discrete probabilities have been considered. Gener- 
alization to continuous distributions is straightforward. The average of a 
function / over phase space T is given by 



(/), = jfd 7 p(7)/(7)- (74) 

The analogue of © for the equilibrium density p( 7 ) is 

p( 1 )=exp K (G((3)-(3H( 1 )). (75) 

The probability-dependent generalization of (J73J) is 

r rp(j) 

«/»,= / d 7 / dxf(j,x). (76) 



r Jo 



Also quantum statistical mechanics has not been considered. The ansatz for 
the equilibrium density matrix p, given a Hamiltonian operator H, is now 

p = ex VK (G(P)-pH). (77) 

It is not obvious how to introduce probability-dependent variables in this 
case. However, for functions of the density matrix, like information content 
/, one can define 

«/» p = [ dxTrpI(xp). (78) 
J o 

It is then straightforward to reformulate the first part of the paper in a 
quantum context. 

Open systems are most often described using stochastic differential equa- 
tions. A recent effort in this direction is found in [2H]- The connection with 
the present formalism has not yet been studied. 
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Appendix A: Ising model 

This appendix explains how the deformed logarithm is constructed in case 
of the one- dimensional Ising model. 

For convenience, let J* = J/ A. Introduce a function X(x), defined for x 
between and 1/2, by the relation 



A(p_(T)) = /3A. 
The inverse function A _1 (z) is given by 



(79) 



sinh(z/2) 



e~ 2J * z + sinh 2 (z/2) 



(80) 



If J = then an explicit expression for X(x) is feasible. One finds in this 
case X(x) = ln(l — x) — \a(x). 

For < x < 1/2 the deformed logarithm must be of the form 



InJx) 



- In (2 cosh(A(x)/2)) - u(X(x)) - 



X(x) 



(81) 



with u(z) still to be determined. The definition for 1/2 < x < 1 is then given 
by 



ln K (x) = ln K (l - x) + A(l - x). 
The derivative of (|RT|) reads 
d 



(82) 



dx 



\n K (x) 



-X'(x) 



L I|l + tanh^) }+u'(X(x)) . (83) 
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Because X(x) is decreasing and the deformed logarithm must be an increasing 
function one obtains the condition that 



2u'(z) + tanh(z/2) > -1. 



(84) 



Because ln re (x) must also increase for 1/2 < x < 1 one has further that 



The condition that ln K (l) = can be analyzed easily. The result is that 
u(z) must vanish in the limit of large z. The condition that ln K (x) is concave 
in x = 1/2 requires that it'(O) > 0. 

The condition of concavity of ln K (x) is more difficult because X(x) is 
necessarily convex for x close to zero, while A"(l/2) = —16 J*. This wrong 
curvature of A(x) close to x — 1/2 should be compensated by the curvature 
of u(z) for z close to zero. 

The above analysis suggests to choose u as follows 



The factor J* has been included to ensure that the deformed logarithm co- 
incides with the natural logarithm in case J* = 0. The above conditions on 
u{z) are satisfied provided J* is less than Jj^ ax ~ 3.83. 

It is rather hard to check concavity of the deformed logarithm in an 
analytic manner. However, by numerical evaluation one can convince oneself 
that the function is concave as long as J* is not too large. 

Appendix B: Thermodynamic stability 

In this appendix is shown that the probability distributions considered in the 
present paper are automatically stable under thermodynamic perturbations. 
It is a tradition (see |T7j) to require also that pressure is a decreasing function 
of volume. However, here the notion of volume dependence is absent. 

Consider the problem of maximizing entropy {{I)) r under the constraints 
that average energy (H) and average magnetization (M) have given values 
U, respectively B: 



2u\z) +tanh(z/2) < 1. 



(85) 



u{z) = l -J\l + z)e-*. 



(86) 



(H) = J2 r k H k = U 



k 
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(M) 



J2 r kM k = B. 



(87) 



Introduce Lagrange multipliers a, (3, and 7, and minimize the expression 
P r kH k + 7 TkMk ~ ^ r ~ a Tk - 

k k k 

Variation w.r.t. r k gives the condition 

= pH k + 1 M k + \n K {r k ) + F K {0)-a. (89) 

Hence the probability distribution maximizing ({I)) r must be of the form 

ln K (r fc ) = -F K (0) + a -(3H k - 7 M fc . (90) 

The parameter a must be such that ^ k r k = 1 and will be considered as a 
function of f3 and 7. The value of ((I)) r at equilibrium is denoted S(P,j). 
Introduce escort probabilities P k (/3,j) by 



1 



z{P,i) 



exp' K (-F K (0) + a-(3H k - jM k ). 



(91) 



where 



Z(P,l) = ]Texp' K (-F K (0) + a -pH h - 7 M fc ). 



(92) 



Then, using (|9~Uj). one obtains 



<9f/ 


k 


da 
_dp 


— Hk 


at/ 

c?7 


k 


da 
d^j 


-M k 


a/? 


= Z(p n )J2PkM k 

k 


da 
_d(3 


— Hk 


c?7 


= Z(p, 7 )J2PkM k 

k. 


da 
c?7 


-M k 



(93) 
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Now, from = 1 follows 



da 

dj3 
da 



where we introduced the notation that (X)* = J2k PkXk- 
This leads to the result that 



dU 
dU 
dB 



-Zfai) [(H%-(H)l] 

dB 
~df3 



— = -Z(/3,7)[(M 2 ),-(M)^]. 

A short calculation using (J§Uj) and (fT^j) gives 

8S_ dU_ dB_ 

dp ~ ^~dp +1 ~dp 

dS_ p^L + d JL 

dj d'j d'-y 



Introduce the notation 



dU dB dU dB 



Then one has, using (IMI) . 

dS_ 

dU 



and similarly 



1 



d(3 c>7 c>7 df3 



dS dB dS dB 



P 



d(3 d"f c?7 dp 



dS 
~dB 



7- 
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Let us now put everything together. The matrix of second derivatives of 
S as a function of U and B equals 



1(100) 



d 2 S d 2 S \ i / 8B _dU 

dU 2 dBdU | - / dy dj 

d 2 S J - D(8 —\ [ _dB dU 

g(/Vr) f (M 2 ), - (M)l -(HM). + (M), 

D(^, 7 ) V ~<™>, + (M)* (if 2 ), - (F) 2 

From Schwartz's inequality follows that 

| - (HM), + (H), (M%\ 2 < [(M 2 ), - (M) 2 ] [(if 2 ), - (ii) 2 ] . (101) 

This inequality suffices to show that the matrix in the r.h.s. of (jlOOj) has 
positive eigenvalues and that D{j3, 7) > 0. Hence, entropy is a concave 
function of U and B. This shows thermodynamic stability of the probability 
distribution r. 
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